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It is widely believed that asphericity in the explosion is the crucial ingredient 
leading to successful core-collapse (CC) supernovae. However, direct 


observational evidence for the explosion geometry and for the connection 
with the progenitor properties are still missing. Based on the thus-far largest 
late-phase spectroscopic sample of stripped-envelope CC supernovae, 

we demonstrate that about half of the explosions exhibit a substantial 
deviation from sphericity. For these aspherical CC supernovae, the spatial 
distributions of the oxygen-burning ash and the unburnt oxygen, as traced 
by the profiles of [Ca 11] AA7291,7323 and [O 1] AA6300,6363 emissions, 
respectively, appear to be anticorrelated, which can be explained if the 
explosion is bipolar and the oxygen-rich material burnt into two detached 
iron-rich bubbles. Our combined analysis of the explosion geometry and 
the progenitor mass further suggests that the degree of asphericity grows 
with the mass of the carbon-oxygen core, which may be used to guide 
state-of-the-art simulations of CC supernova explosions. 


When the nuclear fuel deep inside a massive star (with a zero-age 
main-sequence mass >8 M,) is exhausted, the core collapses to forma 
neutron star or a black hole. A large amount of gravitational energy is 
released, part of whichis transformed into the kinetic energy that drives the 
expulsion of the rest of the star, leading to a core-collapse (CC) supernova. 

The mechanism behind a CC supernova is one of the main 
long-standing challenges in astrophysics. Studying the geometry of 
the explosion and the relation with the mass of the progenitor may shed 
light on this unresolved issue. If the released energy is concentrated 
in a specific direction, the explosion has a bipolar configuration’. 
During this process, the oxygen-rich materials are compressed to the 
equatorial plane, whereas the explosive burning takes place along the 
axis perpendicular to the oxygen-rich plane’*. Bipolar explosions are of 
particular interest due to their potential connection with gamma-ray 
bursts (GRBs). The latter produce a collimated relativistic beam and 
can be observed only when viewed on-axis’. 

The late-phase (nebular) spectroscopy of stripped-envelope (SE) 
supernovae (a subtype of CC supernovae) is a unique tool for studying 


the geometry of the innermost ejecta. An SE supernova is produced by 
amassive star that lost most of its hydrogen envelope before the explo- 
sion®. Several months after the explosion, the ejecta of an SE supernova 
becomes optically thin due to its expansion. As it no longer hidden by 
the massive envelope, the central region is exposed. During this phase, 
the spectrum of an SE supernova is dominated by forbidden lines, in 
particular [O 1] AA6300,6363 and [Ca 11] AA7291,7323 (Fig. 1). Research 
works inthe past focused on [O I] emitted from the unburnt oxygen-rich 
stellar material, the profile of which strongly depends on the viewing 
angle if the explosion is aspherical’*. The double-peaked [O 1], a dis- 
tinct feature of a bipolar explosion, has been found to be common for 
SE supernovae’. This result has been verified with larger samples’. 
However, there is a dilemma when inferring the geometry of the explo- 
sion using only [O 1]: (1) When viewed from the on-axis direction, [O 1] 
will appear to be sharp and single-peaked, making it indistinguishable 
from a spherical explosion. (2) Observing the double-peaked [O 1], 
whichis the smoking gun ofa bipolar explosion, suggests the explosion 
has been viewed on-edge, hiding its possible connection with a GRB’. 
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Fig. 1| Nebular spectroscopy of SE supernovae. The ranges of [O 1] (blue) and 
[Ca 11] (red) are highlighted. From top to bottom: SN 2011dh (SP [O 1] and [Ca 11]), 
SN 2008im (SP [O 1] and horn-like [Ca 11]), SN 2005nb (SP [O 1] and flat-topped 
[Ca 11]) and SN 2004ao (horn-like [O 1] and SP [Ca 11]). The central wavelengths of 
the typical nebular emissions are labelled by the vertical lines. 


(3) Further, an interpretation ofa two-dimensional configuration based 
onaone-dimensional line profile is degenerate” ”?. We, therefore, need 
an improved method to infer the explosion geometry, especially the 
characteristic features when a bipolar explosion is viewed on-axis. 
When a CC supernova explosion occurs, oxygen is burnt into 
iron-peak elements in the direction in which a large amount of the 
energy is released. Outside this region, the material is still oxygen-rich 
and the fractions of the heavy elements are negligible. The explosive 
burning ash and the unburnt material are, therefore, decoupled geo- 
metrically. The relation of the profiles of the emission lines originating 
inthese two regions is, thus, a good probe of the explosion geometry. [O 
I] has been frequently employed to infer the geometry of the unburnt 
oxygen-rich material, whereas [Ca 11], which is strong in most SE super- 
novae, is the main coolant of the oxygen-burning ash, making it anideal 
tool for tracing the configuration of the explosion-made region. 
In this work, we propose combining the analyses of the [O 1] and 
[Ca 11] profiles to infer the geometry of the explosion. Horn-like and 
flat-topped profiles in either of the two lines are considered to be 
signatures of an aspherical explosion. Hereafter, these are collec- 
tively referred to as double-peaked (DP) profiles to distinguish them 
from single-peaked (SP) profiles. We study the nebular spectra of 
80 SE supernovae. Among these objects, 41 are single-peaked [O 1] 
and single-peaked [Ca 11] (OSCaS). For objects considered to have 
experienced a bipolar explosion (that is, with either double-peaked 
[O 1] or double-peaked [Ca 11]), the [O 1] and [Ca 11] profiles appear 
to be anticorrelated. Thus, 27 objects have double-peaked [O 1] and 
single-peaked [Ca 11] (ODCaS) whereas 12 objects have single-peaked 
[O 1] and double-peaked [Ca 11] (OSCaD). Some examples of supernovae 
in different line profile categories are shown in Fig. 1. The anticor- 
relation of the [O 1] and [Ca 11] profiles can be understood as similar 
explosion configurations being viewed from different angles. These 
bipolar explosions are characterized by burning ash distributed along 
a specific axis and are surrounded by an oxygen-rich torus. When 
viewed on-axis, the Doppler shifts of the bipolar calcium elements 
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Fig. 2 | Analysis of the linewidths. a, Relation between the FWHMs of [O 1] and 
[Ca 11]. Objects in the different line profile categories (classified by the profiles 
of [O 1] and [Ca 11]) are labelled by different colours and markers (OSCaS: both 
[O 1] and [Ca 11] are single-peaked; ODCaS: [O 1] is double-peaked and [Ca 11] 

is single-peaked; OSCaD: [O 1] is single-peaked and [Ca 11] is double-peaked). 
Objects with the horn-like [Ca 11] profiles are highlighted by the thick edges. 
Data are presented as mean value with error bars indicating the 16th and 

84th percentiles. b, Statistics of the difference in linewidths. The cumulative 
distributions of the difference in [Ca 11] and [O 1] widths for objects in different 
line profile categories are compared. 


give rise to the horn-like profile of [Ca 11], whereas most of the oxygen 
elements on the torus are moving perpendicularly to the line of sight 
(LOS), producing single-peaked [O 1]. Onthe other hand, when viewed 
froma direction perpendicular to the axis, [O 1] appears to be horn-like 
and [Ca 11] is single-peaked for a similar reason. This simple picture is 
supported further by the data: (1) the line profiles are largely either 
single-peaked or double-peaked and (2) there is no object showing 
both horn-like [O 1] and horn-like [Ca 11]. 

Evidence of asphericity can also be inferred from the widths of 
[O 1] and [Ca 11]. Since the velocity of the material is proportional to 
the radial coordinate under the assumption of homologous expan- 
sion, the width of a line reflects the spatial extension of the emitting 
elements along the LOS”**?°", The full-widths at half-maximum 
(FWHMs) of [O1] and [Ca 11] are compared in Fig. 2a. For the OSCaS and 
ODCaS objects, [O 1] is broader than [Ca 11] in more than 80% of cases. In 
contrast, for OSCaD objects, [Ca 11] is broader than [O 1] in most cases. 
The difference between the [Ca 11] and [O 1] widths is shown in Fig. 2b. 
The OSCaD objects have distinct population from the other twotypes, 
with P< 0.001 from the Anderson-Darling test. The large FWHMs of 
[Ca 11] are not consistent with a spherical explosion in which heavier 
elements are explosively synthesized in inner regions and expelled at 
lower velocities. However, these data can be naturally explained by a 
bipolar explosion in which heavy elements are synthesized along a 
specific axis and ejected at high velocities. If viewed closely from the 
axis, the average LOS velocity of the iron-rich ash may be faster than 
that of the oxygen-rich material, resulting in broader [Ca 11]*. 

The peculiar double-peaked profile of the [Ca 11], along with its 
broadness, strongly suggests a bipolar distribution in the emitting 
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Fig. 3 | Configuration of the idealized axisymmetric model. The O-rich 
region (labelled by blue) and the iron/calcium-rich ash (labelled by red) are 
complementary in their distributions. The profiles given by the model for 
[O 1] (blue) and [Ca 11] (red) strongly depend on the viewing angle, which are 


compared with the observed [O 1] and [Ca 11] profiles (black lines) of SN 2008im, 
SN 2005nb and SN 2004ao. In these panels, the [O 1] and [Ca 11] profiles of SN 
2011dh, an OSCaS object, are plotted to illustrate the difference between the 
spherical and bipolar explosions. 


region. We propose an idealized axisymmetric model to explain the 
profiles ofthe OSCaD and ODCaS objects (Fig. 3). The model is charac- 
terized by two separated zones: the calcium elements are distributed 
ina peanut-like region (Xo = 0 and X,, = 1, where Xis the mass fraction), 
whereas outside this region, the material is oxygen-rich (X, =1 and 
Xca = 0). The physical properties (for example, mass fractions and y-ray 
deposition rate) are assumed to be uniform in each zone for simplicity. 
Noradiation transfer process is included. These simplifications allow us 
to focus on the geometric effect. When viewed from 0 = 15°, 50° or 75° 
from the axis, the model produces single-peaked [O 1]/horn-like [Ca 11], 
single-peaked [O 1]/flat-topped [Ca 11] and horn-like [O 1]/single-peaked 
[Ca 11], respectively, in good agreement with the observed profiles of 
SN 2008im, SN 2005nb and SN 2004ao, which are shown in Fig. 3 as 
examples. It is surprising that sucha simple model involving only two 
parameters (the maximum velocity V,,,, and the viewing angle 6) can 
explain the wide range of both the [O 1] and [Ca 11] profiles simultane- 
ously, as well as their anticorrelation. This result suggests that the 
geometric effect is responsible for the diversity seen in the line profiles. 
Including other ingredients (for example, the density distribution and 
radiation transfer) may improve the agreement with the observations, 
although their roles are probably secondary. 

Statistics is useful for constraining the explosion geometry. The 
occurrence rate of bipolar explosions (Npipotar = Npptcam + Norton) is 0.49 
(=39/80), implying that about half of SE supernovae explode with a 
non-negligible departure from a spherical configuration. Here Nppion 
and Nop ica are the numbers of objects showing double-peaked [O 1] 
and double-peaked [Ca 11], respectively. Within this bipolar explosion 
category, the occurrence rate of double-peaked [Ca 11] is 0.31 (= 12/39), 


which is consistent with the axisymmetric model: at 0 ~ 45°, the [Ca 1] 
profile transforms from double-peaked to single-peaked. If there is no 
bias inthe orientation, the expected occurrence rate of double-peaked 
[Ca 11] in this model is 0.29 + 0.06 for this sample size, assuming that 
the solid angle is uniformly distributed. 

The dependency of the explosion geometry on the progenitor 
mass is crucial for revealing the explosion mechanism of CC super- 
novae. However, observational hints are still missing. In this work, we 
employ the intensity ratio of [O I] to [Ca 11] (hereafter denoted as [O 
1]/[Ca 11]) as a proxy for the CO core mass under the assumption that 
the amount of burning ash is not sensitively dependent onthe CO core 
mass. This relation has been routinely adopted in previous theoretical 
and observational works, with a larger [O 1]/[Ca 11] implying a more 
massive CO core’ >’, To investigate the relation between [O 1]/[Ca 
11] and the asphericity of the explosion, the sample was binned into ten 
groups according to the [O 1]/[Ca 11] sequence. The range of [O 1]/[Ca 
11] for the bins is -0.3 to 0.6 ina logarithmic scale, which roughly corre- 
sponds to a CO core mass of 2.05 to 7.84 M, or azero-age main-sequence 
mass of 12 to 30 M,, based on the scaling relations between these 
quantities’. The relation between [O 1]/[Ca 11] and the bipolar rate 
(Nop tcam + Norton divided by the total number of the objects) in each 
group was investigated, and no correlation was discerned (Spearman 
correlation coefficient p = 0.03), as shown in Fig. 4a. 

The lack of correlation between the bipolar rate and [O 1]/[Ca 11] 
possibly suggests that whether a CC supernovae explodes ina spheri- 
cal or bipolar configuration is not determined by its CO core mass. 
However, the association of [O 1]/[Ca 11] to the CO core mass requires 
further quantification; for example, ifthe explosion energy tends to be 


Nature Astronomy 


Article 


https://doi.org/10.1038/s41550-023-02120-8 


CO core mass (Mo) 


2 3 4 5 6 


7 8 
1.0 1 f 1 1 f L 


Average = 0.49 


Fraction 


— Full sample 
DP [O01] 
= DP [Ca] 


0.8 + 


0.6 5 Z 


8. ans = 60°-65° 


trans 


0.4 + 


Nop [Ca „/NBipolar 


area @,,,,=40°-55° | 


trans 


0.2 4 


6 „n, = 10°-20° 


trans 


o A T T T T T 
-0.6 -0.4 -0.2 o 0.2 0.4 0.6 0.8 
log [O 1]/[Ca 11] 


Fig. 4| The relation between the explosion geometry and the CO core mass of 
the progenitor. a, Occurrence rate of bipolar supernovae as a function of [O 1]/ 
[Ca 11]. The bipolar rate of the full sample (0.49) is indicated by the dashed line. 

b, Cumulative distribution of [O 1]/[Ca 11] for different line profile categories. The 
ratio of the full sample (black), double-peaked [O 1] (pink) and double-peaked [Ca 11] 
(light blue) are compared. c, Occurrence rate of double-peaked [Ca 11] as a function 
of [O 1]/[Ca 11]. The insets schematically illustrate a sequence for the explosion 
geometry asa function of the CO core mass, which explains the observed relation. 
The three models have different transition angels 6,,.,,,. When viewed above this 
angle, the [Ca 11] profile transforms from double-peaked to single-peaked. The 
shaded regions mark the expected double-peaked [Ca 11] rates of these models 
assuming a uniform distribution of the solid angle. In a and c, data are presented as 
mean value with error bars indicating the 16th and 84th percentiles. 


larger for amore massive CO core**”’, thus creating a larger amount 
of Ca (ref. 26,27), this relation may be diluted. In any case, the statisti- 
cal investigation presented here indicates that acommon explosion 
mechanism is behind CC supernovae across a range of CO core masses, 
at least for the sample of canonical SE supernovae studied in the current 
work. Theasphericity of the energy release must be key to producing CC 
supernovae, with about half of the sample having a detectable bipolar 
distribution, as seen for either [O 1] or [Ca 11]. 

A similar analysis was applied separately to the bipolar sample. 
The cumulative distributions of [O 1]/[Ca 11] of the OSCaD and ODCaS 
objects are compared with the full sample in Fig. 4b. Surprisingly, the 
objects with double-peaked [Ca 11] have a larger [O 1]/[Ca 11] ratio than 
those with double-peaked [O I] (P< 0.05 from the Anderson-Darling 
test). Inthe low [O 1]/[Ca 11] range, the presence of double-peaked [O1] 
suggests that the explosion is bipolar. The iron-rich bubbles that char- 
acterize the geometry of the burning ash are not distinctly separated 
into two hemispheres, as inferred from the absence of double-peaked 
[Ca 11]. By binning the bipolar supernovae into five groups according to 
the [O 1]/[Ca 11] sequence, we find a trend in which the occurrence rate 


of the double-peaked [Ca 11] in each group increases with the average 
[O 1]/[Ca 11] (p = 0.87 and P= 0.053; Fig. 4c), although the statistics is 
limited by the relatively small number in each bin (n = 8) and should 
be considered as suggestive”. This tendency can be explained if the 
transition angle 6,,,,,,increases with [O 1]/[Ca 11]. Above this angle, the 
[Ca 11] profile transforms from double-peaked to single-peaked. This 
is consistent with a scenario in which the explosion becomes more 
aspherical and the burning ash becomes more detached and more 
collimated as the CO core grows, which is schematically demonstrated 
in Fig. 4c. 

The correlation of the occurrence rate of double-peaked [Ca 11] 
(Fig. 4c) and the probable non-correlation of the occurrence rate of 
the aspherical explosions (Fig. 4a) to the CO core mass (or the explo- 
sion energy, due to its correlation with the kinetic energy Ex of the 
ejecta’) must be satisfied in the explosion mechanism(s). The stand- 
ard delayed-neutrino explosion mechanism involves various insta- 
bilities’, some of which may introduce a bipolar configuration (for 
example, standing accretion shock instabilities”). Recent simulations 
further show that either rotation or a magnetic field, or both, may have 
an important role in assisting a neutrino-driven explosion’. These 
effects may naturally introduce a specific direction in the explosion 
dynamics. In a massive core, the neutrino-driven mechanism alone 
would be insufficient to overcome the binding energy”, and then the 
rotation and magnetic field may become important even energetically, 
so that there is a magneto-rotational explosion™. 


Methods 

Sample description 

The SE supernova sample in this work is the same as that used in 
ref. 13, which was mainly compiled from the literature*”*®. The analy- 
sis was restricted to nebular spectra acquired not less than 100 days 
after the discovery of the supernova. In this sample, 80 objects have 
spectra with a high enough signal-to-noise ratio to allow an investiga- 
tion of both the [O 1] AA6300,6363 and [Ca 11] AA7291,7323 profiles. 
These were selected for further analysis. Supernovae associated with 
a GRB or X-ray flash (XRF, a low-energy analogue of a GRB), namely 
SN 1998bw/GRB 980425 and SN 2006aj/XRF 060218, were excluded 
from the analysis to avoid introducing any bias in the orientation” *°. 
SN 2012ap, which is suggested to be a GRB-associated supernova but 
viewed off-axis, was also excluded”. The peculiar SN 2005bf, whose 
explosion mechanism is probably different from that of a normal SE 
supernovae, was not included either’**. The SE supernovae studied in 
this work are listed in Supplementary Table 1. 


Linewidth measurement and error estimation 

Before measuring the linewidths of [O 1] and [Ca 11], the nebular spec- 
trum was first corrected for the redshift, which was derived from the 
central wavelengths of the narrow emission lines from the explosion 
site or the redshift of the host galaxy**. The next step was to remove 
the underlying continuum following the same method as used in pre- 
vious works??>*, The line connecting the minima at both sides of the 
emission was defined to be the local continuum and then excised, as 
illustrated in Supplementary Fig. 1. 

To derive the isolated [O 1], we needed to subtract the broad 
Ha-like structure at the red wing of [O I]. This emission was assumed 
to be symmetric with respect to 6,563 A and was constructed by reflect- 
ing the red side to the blue side. A similar procedure was applied to 
derive the [Ca 11] profiles, which were complicated by the many emis- 
sions at similar wavelength ranges emitted by the iron-peak elements. 
However, most of these lines were usually very weak compared with 
[Ca 11] and were not taken into consideration to avoid complications. 
We first subtracted the He1A7065 line by assuming that it was symmet- 
ric with respect to the central wavelength. The profile was constructed 
by reflecting the blue side to the red side. After He I was subtracted, 
the profile of [Fe 1] A7155 was constructed with the same method. 
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The separation of the central wavelengths of [Fe 1] and [Ca 11] was 
~5,700 km s”, whereas the half-width at half-maximum of [Ca 11] was 
always smaller than -5,000 km s” (Fig. 2). The blue side of [Fe 1] was, 
therefore, hardly contaminated by [Ca 11]. It was difficult to remove 
[Ni 11] A7377 from the [Ca 11] complex because its separation from the 
central wavelength of [Ca 11] was small (~3,500 kms”). We assumed 
that [Ni 11] and [Fe 1] have the same profile and that the ratio of the 
intensities L,,7,/L,,,, was 0.6 (for the solar Ni/Fe ratio and temperature 
of the emission region, which was ~4,000 K, see refs. 46,47). The profile 
for [Ca r1] was derived by subtracting the He, [Fe 1] and [Ni 11] profiles 
from the complex. Other contaminating lines were usually very weak 
compared with [Ca 11], and they were not taken into consideration to 
avoid complications. The FWHM of the isolated [O 1] and [Ca 11] were 
then measured and corrected for instrumental broadening. Some 
examples of FWHM measurements are shown in Supplementary Fig. 1. 

The uncertainty in the FWHM was estimated by employing a Monte 
Carlo method. An original spectrum was smoothed by convolving it 
with a boxcar filter. The smoothed spectrum was subtracted from the 
original one, then the standard deviation of the residual flux in the 
wavelength range 6,000 to 8,000 A was estimated as the noise level 
of the spectrum. To generate a simulated spectrum, Gaussian noise 
was added to the spectrum using the estimated noise level. The end 
points of the local continuum, the central wavelengths of the Ha-like 
structure, He 1, [Fe 1] and [Ni 11] were allowed to vary by +600 kms ‘to 
account for the typical resolution of nebular spectroscopy, assuming 
flat distributions. The FWHM was then calculated following the method 
described above. This procedure was repeated 10‘ times, and the 16th 
and 84th percentiles of these measurements were estimated as the 
lower and upper limits of the FWHM. 


Line morphology 

The profiles of [O 1] and [Ca 11], being diverse in morphology, were 
employed to characterize the geometry of the emitting regions. 
Although the [O 1] profile is typically categorized by multi-Gaussian 
fitting, this method has several limitations: (1) The results of the fitting 
can be degenerate, leading to an ambiguous interpretation of the pro- 
file. (2) The blending of emissions from iron-peak elements complicates 
the extension of this approach to [Ca 11]. In this work, we employ the 
convolution method**. The line profiles of [O 1] in SN 2008im and [Ca 
11] in SN 2011dh were employed as templates for single-peaked [O 1] 
and [Ca 11], respectively. These templates were convolved with Gauss- 
ian kernels by freely adjusting the centres and standard deviations 
until their widths matched those of the other objects to ensure proper 
comparison. Profiles that resembled the convolved template were 
classified as single-peaked. Alternatively, profiles that deviated from 
the convolved template but exhibited troughs at -0 km sand roughly 
symmetrical local peaks with respect to 0 km s” were categorized as 
horn-like. Inthe absence of these features, the profiles were classified 
as flat-topped. Note that flat-topped and horn-like profiles could not 
sometimes be unambiguously classified, especially when the troughs 
were shallow. However, the exact distinction between these two profiles 
was not essential to the analysis in this work. Therefore, the horn-like 
and flat-topped profiles were collectively referred to as double-peaked 
profiles. We neglected the substructures that occasionally emerge on 
the broad bases of the emissions and focused on the bulk properties. 
In practice, we found that this evaluation led to unambiguous classi- 
fications between single-peaked and double-peaked profiles. Several 
examples are given in Supplementary Fig. 2. 


Axisymmetric model 

Rather than working on multi-dimensional radiative transfer with 
a detailed treatment of the hydrodynamics and nucleosynthesis, 
we attempted to develop a simplified axisymmetric model that can 
reproduce the unique double-peaked profiles of [O 1] and [Ca 11] and 
their anticorrelation. The axisymmetric model is characterized by 


two separated zones: (1) The explosive burning ash is described by 
two symmetric (with respect to the plane perpendicular to the axis 
of symmetry), off-centre ellipsoids. The velocity of the edge of each 
ellipsoid is Vray the maximum velocity of the model. In this zone, the 
oxygen elements synthesized during stellar evolution are burnt into 
heavier elements, with X, = O0 and X,, = 1 (as we are interested in only cal- 
cium, other iron-peak elements are not included). (2) The O-rich zone 
is enclosed in a sphere with radius Vo = 0.8V,,,,. We used Vo = 0.8 Vnax 
as this value can well produce the widths of both [O 1] and [Ca 11] for 
most supernovae in the sample except for two objects, SN 2006T and 
SN 2012fh. The overlap of this sphere with the iron-rich ellipsoids 
was excised. This zone was filled with the unburnt material where the 
fractions of the iron-peak elements were negligible (Xo = Land Xoca = 0). 
The mass fractions and the density were assumed to be homogeneous. 
Throughout the model construction, we used the velocity as a proxy 
of spatial coordinate. During the nebular phase, the ejecta dynamics 
is well characterized by homologous expansion, such that a grid point 
moves in the radial direction with a velocity proportional to its radial 
coordinate. 

The above model was mapped onto a three-dimensional cube, 
which was divided into 500 x 500 x 500 cells using the Cartesian coor- 
dinates, with each side bounded by -Vna and +V,,,,- No radiation trans- 
fer was introduced. The flux from each grid point was assumed to be 
proportional to its volume. Homologous expansion suggests that the 
infinitesimal volume dr is proportional to its ‘volume’ dV’ in velocity 
space. The specific flux of [O 1] (or [Ca 11]) Fd Vat velocity Vwas calcu- 
lated by summing dV of the O-rich (or iron-rich) grid points with the 
LOS velocities within Vto V+ dV. [O1] is, by nature, a doublet, consisting 
of two transitions emitting at 6,300 and 6,363 A, respectively. These 
two components, separated by -3,000 kms“, have intensity ratio 3:1 
inthe optically thin limit. For the [Ca 11] doublet, the two components, 
centred at 7,291 and 7,323 A, are separated by -1,300 km s”, and their 
intensity ratio was fixed to be 1:1in the optically thin limit. The model 
[O 1] and [Ca 11] profiles, with V,,,,=5,000 km s ‘and viewing angle 0 
varying from 0° (on-axis) to 90° (on-edge), are shown in Supplemen- 
tary Fig. 3. 

The axisymmetric model predicts the unique features of both 
[O 1] and [Ca 11] profiles and their anticorrelation. The [O 1] and [Ca 11] 
profiles are determined by 0 with the opposite tendency. When viewed 
from @ = 0° (on-axis), the model [Ca 11] shows a horn-like profile and [O 1] 
is single-peaked. When viewed from 0 = 90° (on-edge), the [Ca 11] profile 
becomes single-peaked and [O 1] transforms to the horn-like profile. 


Application of the model to the observations 
In this section, the model profiles were applied to supernovae with 
either double-peaked [Ca 11] or [O 1] (except for SN 2015ah, which 
will be discussed later). We freely varied V,,,,, and 0 to find parameter 
values that corresponded to the observations. The procedure was 
as follows. For each object, we focused on the line that exhibited the 
double-peaked feature, that is, [O 1] for the ODCaS category and [Ca 11] 
for the OSCaD category. The viewing angle of the model was varied until 
the depth of the trough in the corresponding line matched the obser- 
vation. Flat-topped profiles were an extreme case in which the trough 
depth was 0. After 6 was fixed, we adjusted V,,, until the widths of the 
observed double-peaked profile and the corresponding model profile 
matched. Any line that appeared to be single-peaked was not directly 
involved in this procedure but served to assess the alignment between 
the model and the observations. In two cases, SN 2006T and SN 2012ah, 
the third parameter, thatis, the maximum velocity of the O-rich region 
Vo, was required, although in most cases fixing ro (= Vo/ Vmax) = 0.8 gave 
reasonable matches. The application of the model to the full sample 
of supernovae that are considered to be bipolar explosions is shown 
in Supplementary Fig. 4. 

The [O 1] and [Ca 11] profiles from the simplified model indeed 
captured the main features of [O 1] and [Ca 11], although they failed to 
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produce some substructures, which may require the introduction of 
clumpy structures (for example, a moving O-rich blob as exemplified 
by SN 1996aq; ref. 11) or a centrally concentrated density profile (SN 
1990aj, SN 2002ap and SN 2006ck). For the horn-like [O 1] profiles, the 
model predicts that the intensity of the peak at the redder wavelengthis 
larger than the bluer one because of the contribution from the second- 
ary component of the [O 1] doublet centred at 6,363 A (or 3,000 kms‘). 
However, some objects had flat-topped profiles (for example, SN1990U 
and SN 2005aj) or a stronger blue peak (SN 2000ew, SN 2004 gt and 
SN iPTF13bvn). This could be explained by the effect of the residual 
opacity: a photon escaping from the far side of the ejecta experiences 
twice the scattering or absorption (see, for example, ref. 20), which 
will effectively reduce the fluxes at redder wavelengths. If the fluxes 
from the grid points with LOS velocity >O were manually reduced by 
10-40%, the model indeed produced [O 1] witha flat-topped profile or 
ahorn-like profile with a stronger blue peak. 

In conclusion, the axisymmetric model can produce the [O 1] and 
[Ca 11] profiles of the sample of SE supernovae with signatures of a 
bipolar explosion used in this work, as well as their anticorrelation. 
Including other ingredients (clumpy structure, centrally concentrated 
density profile or non-axisymmetry) and radiative transfer effects 
may better capture the substructures of the profiles, but the bulk 
properties, that is, a single peak versus a double peak, should not be 
substantially affected. 


Connection with GRBs 

Although we excluded supernovae associated with a GRB from our 
analysis in the main text, the following events are closely examined in 
this section. It has been suggested that SN 1998bw, associated with GRB 
980425, was viewed on-axis, whereas SN 2012ap, also associated witha 
GRB, was viewed off-axis*”**“"”. Existing knowledge on the orientation 
of these two objects allowed us to employ them as test cases to evaluate 
the predictions made by the model developed in this work. 

The [O 1] and [Ca 11] profiles of SN1998bw and SN 2012ap are shown 
in Supplementary Fig. 5. The spectra were observed on 26 November 
1998 (220 days after the explosion”’) and 16 November 2012 (272 days 
after the explosion“), respectively. The [O 1] profile of SN 1998bw 
appears to be single-peaked, whereas that of SN 2012ap has atrough at 
-0 km s”, andis, therefore, classified as double-peaked. Although the 
[O 1] profiles of both objects are aligned with the model predictions, 
their [Ca 11] profiles exhibit peculiar characteristics. For SN 1998bw, 
[Ca 11] appears to be flat-topped witha trough centred at -2,500 km s”. 
The overall shape resembles the [Ca 11] profile of SN 2002ap, a proto- 
type of a flat-topped profile, except for the large blueshift. Ifthe trough 
is viewed as the centre of [Ca 11], the explosion of SN 1998bw must be 
highly aspherical, which cannot be explained by the axisymmetric 
model in this work. In any case, the trough of the double-peaked [Ca 11] 
was predicted to be at -O km s”. The lack of a [Ca 11] flux at the redder 
wing suggests that the burning ash is either highly opaque, thus trap- 
ping the photons from the rear side, or is unipolar*. Additionally, the 
[Ca 11] complex of SN 1998bw was also found to be blended by strong 
[Fe 1] A7155 at its blue wing, which is typically very weak in normal SE 
supernovae. These factors complicate the backward inference of the 
geometry of the explosion-made region from the [Ca 11] profile. Regard- 
ing SN 2012ap, its [Ca 11] is single-peaked but much broader than [O1]. 
A possible explanation for this broadening is contamination by the 
emissions from the iron-peak elements. However, the [Ca 11] of this 
object is too broad for the de-blending procedure employed inthis work. 

To illustrate these points, the model profiles were compared with 
SN 1998bw and SN 2012ap. We fixed the viewing angle to 0 = 0° (SN 
1998bw) and 0 = 75° (SN 2012ap) and then adjusted V,,,,,, until the widths 
of the observed and model [O 1] matched. Although the [O 1] profiles of 
both objects were well reproduced, the model [Ca 11] failed to capture 
the main features of the observed [Ca 11], as illustrated in Supplemen- 
tary Fig. 5. For SN 2012ap, if we assumed that Fe and Ni have the same 


spatial distribution as the Ca and then adjusted the intensity ratios of 
the lines (L7155/L.729, ~ 0.56, L7377/L. 729, ~ 0.12 and Ly4y9/L. 79, ~ 0.31), the 
observed [Ca 11] complex can be reproduced to some extent. Note that 
the model profiles are blueshifted by 1,000 km s” to match the peak 
of the observed profile. However, introducing these emissions cannot 
explain the peculiar [Ca 11] profile of SN 1998bw. In any case, the model 
profiles are too broad. 

The failure in reproducing the [Ca 11] profile of SN 1998bw sug- 
gests that the axisymmetric model may be overly simplified to explain 
the explosion geometry of GRB-associated supernovae. Further, the 
unusually abundance of iron-peak elements seen in both objects indi- 
cates that there were notable differences in explosive nucleosynthesis 
compared to normal SE supernovae. These observations, coupled 
with the absence of a GRB association for normal SE supernovae even 
when viewed on-axis, suggest that a bipolar explosion alone may not 
be sufficient to generate a GRB. More extreme conditions are probably 
required. 


Peculiar object: SN 2015ah 
SN 2015ah is classified as an OSCaD object, with the [Ca 11] complex 
having a horn-like profile with a trough at -4,000 km s7, like the cen- 
tral wavelength of [Ni 11] A7377 (~3,600 km s”). There is also a shallow 
trough at the central wavelength of [Fe 11] A7155 (- -5,600 km s”). 
These features suggest that this complex is probably not dominated 
by [Ca 11], unlike the rest of the sample. Instead, [Ni 11] may be the 
main contributor. Both the [Fe 1] and [Ni 11] lines are emitted from the 
burning ash. Assuming that they share the same spatial distribution 
as the [Ca I1]-emitting region, by adjusting the intensity ratios of the 
lines (L737-/L 79, = 3.7 and L7455/L.799, ~= 1.1), the model indeed reproduces 
the profile of the [Ca 11] complex, as shown in Supplementary Fig. 6. 
The strong [Ni 1] (L7377/L7155 = 3.4) suggests that this object exhibits 
substantial overproduction of the Ni/Fe ratio*®*”. To our knowledge, 
this is the third CC supernova that has been reported to have sucha 
super-solar Ni/Fe ratio. The other two cases are SN 2012ec (ref. 47) and 
SN 2006aj (ref. 50). Although this sample is very small, it seems that 
an excess of the Ni/Fe ratio can occur for all types of CC supernova, 
irrespective of the presence or absence of the outer envelope. Further, 
two out of the three cases, SN 2006aj and SN 2015ah, are likely associ- 
ated with non-spherical explosions, as an asymmetry in the explosion 
can effectively eject the deep-lying burning ash, leading to the change 
inthe Ni/Fe ratio’. A systematic investigation of the relations between 
the Ni/Fe ratio and the properties of the explosion will be important 
to reveal the origin of the iron-peak elements and provide constraints 
onthe explosion mechanism of CC supernovae. 


Data availability 

Most of the spectra are available from WiseRep (https://www.wiserep. 
org/) and Supernova Database of Berkeley (http://heracles.astro.berke- 
ley.edu/sndb/). The data that support the plots within this paper and 
other findings of this study are available from the corresponding author 
upon reasonable request. 


Code availability 

Astropy, Matplotlib, Numpy and Scipy are available from the Python 
Package Index (PyPI) (https://pypi.org/). Upon request, the first author 
will provide the Python code used to generate the model line profiles. 
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